r]- & r/'-mesic nuclei and Ua{^) anomaly at finite density 
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We discuss theoretically the possibility of observing the bound states of the ri and 77'(958) mesons 
in nuclei. We apply the NJL model to study the r) and 77' meson properties at finite density and 
calculate the formation cross sections of the rj and 77' bound states with the Green function method 
for (7,p) reaction. We also discuss the experimental feasibility at photon facilities like SPring-8. The 
contributions due to the u) meson production are also included to obtain the realistic (7,p) spectra. 
We conclude that we can expect to observe resonance peaks in (7,p) spectra for the formation of 
meson bound states and we can deduce new information on 77 and rj' properties at finite density. 
These observations are believed to be essential to know the possible mass shift of r)' and deduce new 
information on the effective restoration of the ?7a(1) anomaly in the nuclear medium. 

PACS numbers: 25.20.-x, 14.40.-n, 36.10.Gv 



I. INTRODUCTION 

Understanding the low-lying liadron structures from 
the viewpoint of the quark and giuon degrees of freedom 
is one of the most challenging problems in the quantum 
chromodynamics (QCD) because of its nonperturbative 
nature in the low-energy regime. 

The lightest excitation on the QCD vacuum is the pion 
which is considered as a quark and an antiquark bound 
state in the pseudoscalar channel. Its mass is off-scale 
light compared with other hadrons. This can be under- 
stood by recognizing that the chiral symmetry is spon- 
taneously broken in the QCD vacuum and the pion is 
the Nambu-Goldstone (NG) boson associated with the 
dynamical chiral symmetry breaking (DCSB). It is be- 
lieved to be responsible for a large part of the constituent 
quark masses, which are introduced in the many con- 
stituent quark models. The effects of the explicit break- 
ing of the chiral symmetry on the pion properties have 
been systematically studied using the effective lagrangian 
composed of the pion field. This approach is called the 
chiral perturbation theory (ChPT) The success of 
ChPT approach supports the importance of the DCSB 
in low-energy QCD. 

When we start looking at the strange quark sector, 
we encounter another problem. The ninth heavier pseu- 
doscalar meson is rj' and its mass is much heavier than 
other octet pseudoscalar mesons. Weinberg showed that 
the mass of r]' should be less than J^rriT^ if J7^(l) sym- 
metry were not explicitly broken 0]- Thus the Ua{^) 
symmetry must be broken. The key step to solve this 
C/^(l) problem was to realize that C/4(l) symmetry was 
explicitly broken by the quantum anomaly. In the fol- 
lowing year, 't Hooft pointed out the relation between 
anomaly and topological gluon configurations of 
QCD, i.e., instantons and showed that the interaction 
of light quarks and instantons breaks the Ua{^) symme- 
try He also showed that such an interaction can be 
represented by a local 2Nf quark vertex, which is anti- 
symmetric under flavor exchanges, in the dilute instanton 



gas approximation. 

The effects of the Ua{^) anomaly on the low-energy 
QCD have been extensively studied in the 1/iVc expan- 
sion approach 0. In the Nc ^ 00 limit, the C/yi(l) 
anomaly is turned off and then the 77 and 77' mesons 
become the ideal mixing states. The higher order ef- 
fects of the 1 /Nc expansion give rise to the flavor mixing 
between the 77 and 77' mesons and push up the 77' mass 
0, 0. 01 ■ They were further discussed in the context of 
the ChPT 1^ and the reasonable description of the nonet 
pseudoscalar mesons was obtained. 

The J7yi(l) breaking 2Nf quark determinant interac- 
tion was introduced to the low-energy effective quark 
models of QCD. The low-l ying meson properties have 
been studied [IS [H El El El El El using the three- 
flavor version of the Nambu-Jona-Lasinio (NJL) model 
[13 with the Kobayashi-Maskawa-'t Hooft (KMT) de- 
terminant interaction [1 0| . It is further argued that 
the J7A(l)-breaking interaction gives rise to the s pin-spin 
forces, which are important for light baryons 19l l2Cll21| . 

It was observed by Witten that the instanton contri- 
bution scales as exp(— iVc) which seems to contradict the 
1/Nc expansion approach, m^, = 0{1/Nc) A few 

years ago, Schafer showed that the large Nc behavior 
of the instanton approach is consistent with the 1/Nc 
expansion approach if the instanton ensemble is stabi- 
lized by a classical repulsive core The use of the 
[/^(l) breaking quark determinant interaction with the 
low-energy effective quark model of QCD in the mean 
field approximation seems to mimic the dilute instanton 
liquid picture, and therefore, we consider this approach 
is not inconsistent with the 1/Nc expansion approach. 

The dynamics of instantons in the multi-instanton vac- 
uum has been studied by many authors, either in the 
models or in the lattice QCD approach, and the widely 
accepted picture is that the QCD vacuum consists of 
small instantons of the size about 1/3 fm with the den- 
sity of 1 instanton (or anti-instanton) per fm"*^ In 
the instanton liquid model, the instanton plays a crucial 
role in understanding not only the Ua{1) anomaly but 
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also the spontaneous breaking of chiral symmetry itself. 

In this paper, we consider the rj and rj' meson bound 
states in finite nuclei, i.e., 77- and ry'-mesic nuclei in order 
to investigate the change of the QCD vacuum structure 
at finite density through the 77 and 77' meson properties 
in the nuclear medium. Since the vacuum properties of 
1] and 77' are believed to have close relation to the QCD 
vacuum structure as described above, we can expect to 
have new information and much deeper insights on the 
QCD vacuum by knowing the 77 and 77' meson proper- 
ties at finite density, where the partial restoration of the 
DCSB is expected. 

The 7;-mesic nucleus was studied by Haider and Liu [25l | 
and Chiang, Oset and Liu ,2^]. As for the formation re- 
action, the attempt to find the bound states by (7r"'",p) 
reaction led to a negative result |23|. Recently, the 77 
meson and light nuclei systems are studied experimen- 
tally and also theoretically 529|. The experimental 
study in Ref. [2^ claims the existence of the 7;-bound 
state in the light nucleus. 

Recently, there are several very important develop- 
ments in the research of the spontaneous breaking of 
chiral symmetry and its partial restoration at finite den- 
sity by studying the hadronic systems, such as pionic 
atoms j32j , 77- and 7;'-mesic nucleiJsSJH HE 

El m and w-mesic nuclei El El EflETin both 
of theoretical and experimental aspects. There is also a 
theoretical study on a meson nuclear bound states based 
on SU(2) hnear sigma model |42|. Especially, after a se- 
ries of deeply bound pionic atom experiments [isl , 
K. Suzuki et al. reported the quantitative determina- 
tion of pion decay constant /^r in-medium l31l from the 
deeply bound pionic states in Sn isotopes |45j and stim- 
ulated many active researches of the partial restoration 
of chiral symmetry at finite density Hfl HI E Elii ■ 

However, as for the behavior of the Ua{^) anomaly 
in the nuclear medium, the present exploratory level is 
rather poor. Although some theoretical results have been 
reported, there exists no experimental information on 
the possible effective restoration of the Ua{^) anomaly 
at finite density. In the context of the instanton dy- 
namics, the mass of the 77' at a nonzero baryon chem- 
ical potential was studied in the two color QCD [i^ . 
As for the case of the finite temperature, the tempera- 
ture dependence of the instanton density was calculated 

m "50", 17, m mill. In the low-energy effective model 

of QCD approach, T. Kunihiro studied the effects of the 
L''a(1) anomaly on 77' properties at finite temperature us- 
ing the NJL model with the KMT term and showed the 
possible character changes of 7/' at T 7^ [s^. Theo- 
retical predictions by other authors also supported the 
possible change of the 77' properties at finite density as 
well as at finite temperature ^5S , .57j . Experimental fea- 
sibilities to observe the r;'-mesic nuclei were also studied 
in this context in our previous work in Ref. |37|. where 
the in-medium rj' properties were estimated quite roughly 
based on the results reported in Refs. [s^ [s^ In 
Ref. [3^ further theoretial studies have been performed 



for 77- and 77'-nucleus systems with paying attention to 
the QCD axial C/(l) problem. 

In this study, we consider the rj and t]' meson properties 
in atomic nuclei and the structures of the 77- and 77'-mesic 
nuclei using the SU(3) / NJL model in order to get much 
deeper insights on the 77' mass shifts and their relation 
to Ua{^) anomaly effects and to get more quantitative 
results than our previous work 37]. We also propose 
the formation reaction of the 77- and 77'-mesic nuclei and 
discuss the possibility to produce the rj- and 7/'-nucleus 
bound states. The rj and 77' properties in the medium 
should provide us important information on the QCD 
vacuum structure, especially on the effective restoration 
of the Ua{^) symmetry in the nuclear medium. 

This paper is organized as follows. In SecE] we will 
obtain the quark masses and the quark condensates in 
finite density using NJL model, and calculate the meson 
masses in-medium. In Sec lIIII we will discuss the forma- 
tions of mesic nuclei using the Green function method, 
and show the numerical results of the missing mass spec- 
tra of the (7,p) reaction. Finally, Sec llVI is devoted to 
summary of this paper. 



II. QUARK AND MESON MASSES IN NJL 
MODEL 



We work with the following NJL model lagrangian den- 
sity extended to three-flavor case: 



Co 



— Co + C4 + Cq, 

— tp {idfij^ — rh) Tp . 



9s_ 
2 



^ [(7^A»' + (7AA'^n57^)' 



Ca = .g_D {dot [7/5^(1 - 75)'0j] + /i.e.} 



a=0 



(la) 
(lb) 

(Ic) 

(Id) 



Here the quark field t/j is a column vector in color, fla- 
vor and Dirac spaces and A" (a — . . . 8) is the Gell- 
Mann matrices for the flavor ?7(3). The free Dirac la- 
grangian £0 incorporates the current quark mass matrix 
777 — diag(77i„, 777,^, 771s) which breaks the chiral {7^(3) x 
Ur{?i) invariance explicitly. £4 is a QCD motivated four- 
fermion interaction, which is chiral Ul{^) x Ur{3) in- 
variant. The Kobayashi-Maskawa-'t Hooft determinant 
interaction Cq represents the Ua{^) anomaly. It is a 3 x 3 
determinant with respect to flavor with i,j = u, d, s. 

Quark condensates (qq) and constituent quark masses 
Mq are self-consistently determined by the gap equations 
in the mean field approximation, 

Mu = ruu- 2gs (uu) - 2gD (dd) (ss) , 
Md = rud- 2gs (dd) - 2gD {ss) {uu) , 
Ms = rus- 2gs {ss) - 2gD {uu) (dd) , (2) 
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with 



with 



Tr^"'^) [iSU^ = 0)] 



-2N, 



(2^)4 



Ma 



(27r)3 



(3) 



Here, the uhraviolet cutoff A is introduced to regularize 
the divergent integral and Tr'^'^'^^ means trace in color 



and Dirac spaces. Eq — \Jp\ + Mq is the on-shell energy 

of the quark and rip and rip are Fermi occupation numbers 
of quarks and antiquarks, respectively, defined as, 



npiT,tiq) = J 
np{T, fiq) 



1 



1 



1 + e(^<!+M<j)/5 



(4) 



where T represents the temperature of the system and 
fiq is the quark chemical potential. Up and rip depend on 
the momentum Pq through the energy Eq. For T = 0, 
we can simply write Eqs. Q as 



np{0,IJ,q) = 9{flq-Eq), 

npiO,Hq) = 0, 



(5) 



and by integrating them in momentum space we obtain 
the quark number density pq as 



\[^,l-Mlf/^9{p,q-Mq). 



(6) 



To simulate the symmetric nuclear matter, we consider 
SU(2) symmetric quark matter defined as pu — Pd, Ps = 
0. In this case, the nucleon density p is defined as. 



(7) 



where pu = Pd and Af„ — Md- The NJL model at finite 
temperature and density has been reviewed in |58j| . 

The pseudoscalar channel quark-antiquark scattering 
amplitudes, 

(p3,P4;OUt \pi,P2;m) = {2TT)^6^{p3+P4-pi-p2)Tqq (8) 

are then calculated in the ladder approximation. We as- 
sume that rriu = rud so that the isospin symmetry is 
exact. In the rj and rj' channels, the explicit expression 
of Tqq is 



■'qq — 



u{p3)X^v{p4,) 

u{p3)X°v{pi) 

v{p2)X^u{pi) 
v{p2)X°u{pi) 



Aiq^) B{q^) 



(9) 



A{q') = 
C{q') = 



detD(g2 
2 



detD(g2) 
2 

detD(g2) 



y{2G/"(g2)-G8}, (10a) 
{-2Gr"(g2)-G„}, (10b) 
{2G/«(a_Go}, (10c) 



where q = pi + P2 and G = GqGs ~ GmGm, and 

Go = - ^(2(mu) + (ss))5D, (Ha) 



Gs = -^gs - g(('5s) - Muu))gD 



Gr 



3\/2 



((ss) - {uu))gD ■ 



(lib) 
(11c) 



The quark-antiquark bubble integrals are defined by 

r d^p 



(27r)4 



X Tr(='^'^) [Sf{p)X°i15Sf{p + g)A°i75] , (12a) 



/»(g- 



_(Pp_ 

(2^ 



X Tr(='/'^) [SF{p)X''tj5SF{p + 9)A««75] , (12b) 



d'^p 



(27r)4 

X Tr'^'-^'^^ [SF{p)X°lj5SF{p + q)X''lj5] 



(12c) 



In the calculations of the quark loop integrals, we have 
applied standard techniques of thermal field theory, i.e., 
the imaginary-time formalism and we have introduced 
the same three-momentum cutoff used in the gap equa- 
tions lEI and ©. 

The 2x2 matrix D is given by 



with 



^^'^ > - [ D2iiq') D22{q') 



Dniq') = 2Gol'{q') + 2G„J"\q^) 

D,2iq') - 2Gor"(g') + 2G„/«(g2) 

D2iiq') - 2GsI"\q') + 2GmI%q') 

D22{q^) = 2Gsl'iq') + 2G^r\q') 



(13) 



1, (14a) 
(14b) 
(14c) 

1 . (14d) 



From the pole positions of the scattering amplitude Eq. 
0, the 77-meson mass and the ry'-meson mass m^/ 
are determined. 
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The scattering amplitude Eq. @ can be diagonalized 
by the rotation in the flavor space 



J-qq 



u(p3)A^w(p4) 

u(p3)A°w(p4) 

V{p2)>?u{pi) 

u(P2)A°m(pi) 
u{p:i)yi v{p4) 

v{P2)y^u{pi) 

v{p2)X'^'u{pi) 







B{q^) C{q^) 





D'^'iq^) 



(15) 



with A'' = cos6'A« - sin6'A^ A**' = sin 61 A^ 



cos( 
sin^ 



— smf 
cos 6 



The rotation angle 9 is determined by 

me) 



tan 26* = 



cos ^A*' and 
(17) 

(18) 



Note that 9 therefore depends on q^ . At q^ = iri^, 9 
represents the mixing angle of the A* and A° components 
in the jy-meson state. In the usual effective pseudoscalar 
meson lagrangian approaches, the rj and 77' mesons are 
analyzed using the g^-independent rj — rj' mixing angle. 
Because of the g^-dependence, 9 cannot be interpreted as 
the T} — rj' mixing angle. The origin of the g^-dependence 
is that the rj and 77' mesons have the internal structures. 

We obtain the dynamical quark and meson masses in 
vacuum as compiled in Table ^with the input parameters 
determined in Ref. (53] . In Ref. [s^ , the 4 model param- 
eters, namely, cutoff A, 4-quark coupling constant gs, 6- 
quark coupling constant g^) and the current s-quark mass 
rris have been fixed so as to reproduce the observed values 
of the pion, kaon, 77' masses and the pion decay constant, 
while the current u, d-quark mass has been fixed at 5.5 
MeV, which has been taken from the results of the chiral 
perturbation theory and QCD sum rule approaches. This 
parameter set gives an rj mass of ttt,^ = 514.8 MeV, which 
is about 6% smaller than the observed value = 547.75 
MeV. One is able to fit the mass of 77 instead of the 77'. 
The reason why we have used the present parameter set is 
just we are more interested in the 77'-mesic nucleus than 
the 77-mesic nucleus in this article because we consider 
that the former is more suitable for observing the finite 
density effect of the Ua{^) anomaly. 

In this paper, we investigate the finite density effects 
on the meson mass spectra for the following 3 cases with 
the different strengths gn of the determinant KMT in- 
teraction as. 



(a) gD{p) = go 

(b) go^p) - 

(c) goip) = gDexp[-(p/po)2] 



(19) 



Input parameters 


Calculated results [MeV] 


A = 602.3 [MeV] 


M„,d = 367.6 


u -- 


= 92.4 


gsA^ = 3.67 


M, = 549.5 




= 135.0 


goA^ = -12.36 


{uuY/^ = -241.9 


777,, 


= 514.8 


7n„,d = 5.5 [MeV] 


[ss)'^''-^ = -257.7 


777^/ 


= 957.7 


ms = 140.7 [MeV] 









TABLE L Input parameters determined in Ref. |59ll and cal- 
culated results in vacuum in three-flavor NJL model. 
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FIG. 1: Density dependence of the quark condensates in the 
SU(2) symmetric matter, p„ — pd and ps — 0, where (uu) = 
(dd). Three panels correspond to the cases (a), (b) and (c) 
defined in Eq. 119II . respectively. The nucleon density p is 
defined in Eq. Q and po is the normal nuclear density po = 
0.17 fm-^ 



where go is the vacuum strength of the determinant in- 
teraction as shown in Table The goip) has no density 
dependence for cases (a) and (b). In case (a), the me- 
son vacuum properties are well reproduced as shown in 
Tabled while there are no anomaly effects in case (b). 
For gu — Q case, we use sli ghtl y different parameter set 
as shown in Table I in Ref. W| to reproduce the meson 
masses and the pion decay constant in vacuum without 
anomaly effect. In case (c), we simply assume the density 
dependence of go as this form in order to examine the 
medium effect due to density dependence of go itself on 
the meson mass spectra in finite density. 

Here it may be interesting for our study to notice that 
there are theoretical s^uggestions about possible density 
dependence of [6lll62| \ In Ref. the effective cou- 
pling constant of the instanton-induced interaction is sug- 
gested to have chemical potential dependence for Nf = 2 
systems. For Nf = 3 systems, we can expect to have the 
similar ji dependence, though it is not easy to show ex- 
plicitly. We are interested in studying the effect of such 
density dependence discussed in Ref. j62| on meson mass 
spectra as future works. 

In Fig. n we show the calculated quark condensates 
as functions of density for three types of go (p) defined 
in Eq. (|19|l . In case (b), we have switched off the effect 
of the instanton-induced flavor mixing interaction and 
therefore the s-quark condensate has not changed in the 
SU(2) symmetric matter. The absolute values of the u,d- 
quark condensate decreases significantly faster in case (c) 
than in case (a) when the density goes up. It means that 
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FIG. 2: Density dependence of the meson mass spectra. Three 
panels corresponds to the cases (a), (b) and (c) defined in 
Eq. I19II . respectively. The nucleon density p is defined in 
Eq. ||7|l and po is the normal nuclear density po ~ 0.17 fm~^. 



the contribution of the instanton-induced interaction on 
the u,d-quak condensate is sizable. 

In Fig. 13 we show the calculated density dependence 
of the meson mass spectra for three cases defined in 
Eq. H19|) . In the case of constant go (a), we find that 
the mass of rj' decreases rapidly as a function of the den- 
sity, while the masses of tt and r/ gradually increase. The 
instanton-induced interaction is repulsive for the flavor 
singlet qq channel and the effective coupling strength is 
gD{{uu) + (Jd) -I- (ss))/3. Since the absolute values of 
the quark condensates decrease as the density increase, 
the effective repulsive interaction in the flavor singlet qq 
channel becomes small as the density increases. That is 
the reason why the i]' mass decrease. In Fig. OJb), we 
find that tt and rj are degenerate completely and their 
masses increase gradually as density, and the mass of ij' 
has no density dependence without the ?7yi(l) anomaly 
effects. Without the Ua{^) anomaly effects, the rj' is 
the ideal mixing state, i.e., the pure ss, and therefore no 
density dependence in the SU(2) symmetric matter with 
fig = 0. It should be mentioned that the masses of rj and 
rj' at p = are not reproduced in this case. As shown in 
Fig.EIc), if the go has density dependence as case (c) in 
Eq. (|19|) , we can expect mass reduction of both rj and rj' 
mesons up to around p = po. 

It is interesting to compare our results in Fig. |21 with 
those obtained in Ref. [s^ , where the Quark-Meson Cou- 
pling model (QMC) was used to obtain the mass-shift for 
given mixing angles. The results in Table 1 of Ref. [s^ l 
show that the both masses of rj and rj' decrease in medium 
at normal nuclear density (defined as po = 0.15 fm~^ 
in Ref. ^]) for all cases considered there. These re- 
sults are not consistent to our results of case (a) where 
the 7] mass slightly increases at the normal nuclear den- 
sity. Our results of case (c) are qualitatively consistent 
to those in Ref. [s^ and predict the mass reduction of 
77 and 77' mesons in medium at normal nuclear density, 
however, the sizes of the mass shift are significantly dif- 
ferent in these results. The origin of the discrepancies 
in theoretical calculations are not clear, however, we be- 
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FIG. 3: Density dependence of the quark condensates (left 
panel) and the meson masses (right panel) are shown for the 
SU(2) symmetric matter (thick lines) and the SU(3) symmet- 
ric matter (thin lines). The nucleon density p is defined in 
Eq. Q and po is the normal nuclear density po = 0.17 fm~^. 



lieve that this fact indicates the significant importantce 
of experimental information. 

Here, in order to see the effects of the species of the 
matter to the calculated results, we compare the results 
in the SU(2) symmetric matter with those in the SU(3) 
symmetric matter, pu — pd — Ps, ^''^ Fig. |31 In Fig. |31 
the strength of the determinant interaction is fixed to be 
constant as case (a) in Eq. (fT^ . We get same results as 
those obtained in Ref. [63 in SU(3) case, and find that 
the differences between the calculated results in SU(2) 
and that of SU(3) are less than 10% at normal nuclear 
density pQ. 



III. FORMATION SPECTRA OF MESIC 
NUCLEI 

In previous section, we discuss the possible changes of 
the meson masses due to the medium effects. If such mass 
changes occur in the nuclear medium, we can translate 
these effects into the potential forms. We estimate the 
real parts of the rj- and Ty'-nucleus optical potentials based 
on the calculated their masses at finite density by the 
NJL model. The optical potential U{r) can be written 
as. 



U{r) = V{r)+iW{r), 



(20) 



where V and W indicate the real and the imaginary parts 
of the optical potential, respectively. The mass term in 
the Klein-Gordon equation for the 77 and rj' mesons can 
be rewritten at finite density as; 



ml 



inn?{p) = (mo -I- Am(p))^ 
^ ttiq -|- 2?7ioAm(p), 



(21) 



where mg is the calculated meson mass in vacuum and 
m{p) the meson mass at density p. The mass shift Am(p) 
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is defined as Am(/9) = m{p) — itiq and is Am(0) = 0. 
Thus, we can interpret the Am(p) as the strength of the 
real part of the optical potential, and in this paper we 
write V{r) as 



V{r) = Am(po) 



Po 



(22) 



using the meson mass shifts at normal nuclear density 
Pq. Here p{r) is the nuclear density distribution, which 
is assumed to be of an empirical Woods-Saxon form as 



p{r) = 



Po 



1 + exp 



(23) 



where R = 1.18^43 - 0.48 fm and a = 0.5 fm with the 
nuclear mass number A. 

From the calculated results shown in Fig. [21 we can 
deduce the real part of the rj- and ry'-nucleus optical po- 
tential as defined in Eq. (|22ll . We show the real parts 
of the optical potentials in Fig. 0] for all these cases de- 
fined in Eq. H19(l . We find that the sign and strength 
of the real part of the optical potential strongly depends 
on the anomaly effects at finite p. The potential of rj' is 
attractive for cases (a) and (c), and its strength depends 
on the density dependence of goip)- As we can see in 
Fig. EJb) , no real potential is expected for ij' meson for 
go = Q case. As for ry-nucleus optical potential, even 
the sign of the real potential is changed by the density 
dependence of goip) as shown in Figs.E^a) and (c). For 
no anomaly case (b), "q meson degenerates with pion and 
the real potential is slightly repulsive. Thus, we confirm 
that we can expect to obtain the information on £^4(1) 
anomaly by knowing the t] and ry' meson properties at 
finite p. 

The strength of the imaginary potential of meson and 
nucleus is estimated separately since it is extremely diffi- 
cult to evaluate it using NJL model. For the r] meson, we 
estimate the strength of the imaginary potential based on 
the results obtained in Refs. |3g| assuming the iV*(1535) 
dominance in the rjN channel, and take as; 



W{r) 



-40^[MeV]. 
Po 



(24) 



As for the -q' nucleus optical potential, we can esti- 
mate W{r) from analysis of 7p — > -q'p data fg^l. Since 
they included only A^*(1535) as a baryon resonance in 
the analysis of the rj' formation reaction and determined 
?7'A^A^*(1535) coupling strength, we can easily calculate 
the 77' self-energy in the medium in analogy with the A- 
hole model for the 7r-nucleus system as. 
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P 



2m,y rriri' + Mn 



(+77 -8i)— [MeVl. 

PO 



We consider the reasonable running range is —5 MeV 
- — 20 MeV for the strength of the imaginary part of 



50 



-50 

> -100 

> -150 
-200 
-250 
-300, 



ri (a) 


n (b) 


(c) 






h' 



2 4 
r[fm] 



2 4 
r[fm] 



2 4 
r[fm] 



FIG. 4: Real parts of the r\- and 77'-^^B optical potentials de- 
termined by the meson mass shifts of NJL model in SU(2) 
symmetric matter as defined in Eq. I22II . Three panels cor- 
respond to the cases (a), (b) and (c) defined in Eq. 1191 1. 
respectively. 



the 77' nucleus optical pot ential based on this evaluation 
in Eq. as in Ref. |33- The ry'-nucleus bound states 
were calculated theoretically before in Ref. where the 
widths of the Ty'-mesic nuclear states were not evaluated. 

We should mention here that the evaluation in Eq. H25|) 
provides the repulsive real part, which is opposite to the 
evaluation from the ry' mass shift, since we consider only 
TV* (1535) resonance with the mass Mat. < to^/ -|- A/at 
as an intermediate states. By the (7,p) experiments pro- 
posed in this paper, we can expect to distinguish these 
potentials and to determine the sign and strength of the 
77'-nucleus optical potential. 

In order to observe meson-nucleus bound states exper- 
imentally, we consider missing mass spectroscopy 65], 
which was proved to be a powerful tool for the meson 
bound states formation in the studies of deeply bound 
pionic states [6^ . In this spectroscopy, one observes 
an emitted particle in a final state, and obtains the 
double differential cross section d^a/d^/dE as a func- 
tion of the emitted particle energy. As discussed in 
Refs. li^, E3j EE E^i we think that the (7,p) reaction 
with GeV photon beam is one of the appropriate reac- 
tions for our purpose and it can be performed in existing 
facilities hke SPring-8. We consider the (7,p) reactions 
for the 77- and T^'-mesic nuclei formation. 

We choose the incident photon energy as = 2.7 
GeV, which is the beam energy accessible at SPring-8, 
and choose ^^C as a target nucleus as in Ref. JST^. We use 
the Green function method (7^] to calculate the formation 
cross sections as, 



\ dVldE J 



A{-t,p)ri,ri'tg,{A-l) 



d(T\ 



Lab 



dnj^ 



]SiE), 



(25) In 



where S{E) 

Lab 



(26) 

is the nuclear response function and 



, are the elementary cross sections in the lab- 
oratory frame for the 77 and 77' meson production. The 
elementary cross section for rj' production is estimated to 
be 150 nb/sr using the data of SAPHIR collaboration [tJI 
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FIG. 5: The calculated spectra of the ^^C(7,p)"B® rj' reac- 
tion at — 2.7 GeV are shown as functions of the excited 
energy Eex defined in the text. Eq is the rj' production thresh- 
old energy. The real part of the r;'-nucleus optical potential is 
evaluated for case (a) in Eq. The imaginary part of the 

potential is assumed to be W{r) — —5p{r)/po [MeV] (upper 
panel) and W{r) = -20p(r)/po [MeV] (lower panel). The 
total spectra are shown by the thick solid line, and the dom- 
inant contributions of subcomponents are shown by dotted 
and dashed lines as indicated in the figures. 



and their analysis 73 • As for the i] production process, 
we take the same value for the elementary cross section 

^ (^)p(7p)77 ^ nb/sr based on the data at SPring- 
8 llSl which show the similar production rate for both 
mesons in a test experiment in the preparation stage. 
We sum up all (proton-hole) (g) (meson- particle) configu- 
rations in the final state to get the total cross section in 
Eq. • 

To calculate the response function S{E), we use the 
Green function G{E;f,r') defined as fK^ . 



case (a) 


(fl 

T? A 1 


\ ^' 




case (b) 


n J" 






case (c) 


CO 
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FIG. 6: 

for the ri, 



The calculated spectra of the ^^C(7,p) reactions 



UJ, r) 



mesic nucleus formation at 



2.7 GeV 



are shown as functions of the emitted proton energies in final 
states. The three cases (a), (b) and (c) for go defined in 
Eq. 11911 are considered. The vertical dashed lines indicate 
the production thresholds of the -q, uj and r)' mesons. The 
contributions from the 77, uo and Ty'-mesic nuclei are shown 
by dotted, dash-dotted and dashed lines, respectively. The 
imaginary part of the ry'-nucleus optical potential is W{r) = 
— 20p(r)/po MeV. The oj-nucleus optical potential obtained in 
Ref. E3l is used and is attractive. 



> 
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G{E-r,?)^{p-^\<l,{r) 
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where 0^ is the meson creation operator and \p^^) is 
a proton hole state. The Hamiltonian H contains the 
meson- nucleus optical potential U. We can rewrite 
Eq. H27|) in a simple expression as [t^ , 



G{E- f, ?) ^ YiUmMr')GiiE; r, /) 

1. 771 

Gi{E;r,r') = -2mfcu/(fc, r<)w|+^ (A:, r>). 



— (?!)■'■ (r')lp"^), (27) FIG. 7: The calculated spectra of the ^^C(7,p) reactions 

for the rj, u, 77' mesic nucleus formation at E-y = 2.7 GeV 
are shown as functions of the emitted proton energies in final 
states. The three cases (a), (b) and (c) for go defined in 
Eq. 11911 are considered. The vertical dashed lines indicate 
the production thresholds of the 77, uj and rj' mesons. The 
contributions from the 77, uj and 77'-mesic nuclei are shown 
by dotted, dash-dotted and dashed lines, respectively. The 
imaginary part of the Ty'-nucleus optical potential is W(r) = 
— 20p(r)/po MeV. The oj-nucleus optical potential obtained in 
Ref. )4l|| is used and is repulsive. 



(28) 
(29) 



8 



where ui and are the radial part of the regular and 
outgoing solutions of equation of motion, respectively. 
Using the Green function, the response function can be 
expressed as 

SiE) = --Im f d^rdad^r'da' 

X f{r,a)G{E;r,r')f{r',a). (30) 
Here, we define /(r, a) as 
/(r, a) = x}im ^M) [yrir) ® V-j, (r, a)] 

. 

where Xi s-^d Xf denote the projectile and the ejectile 
distorted waves, respectively, is the proton-hole wave- 
function and ^ is the spin wavefunction introduced to 
count possible spin directions of the proton in the target 
nucleus. The distorted waves Xi Xf are evaluated by 
using the eikonal approximation as in Ref. |37| |. 

In Fig. El we show the calculated spectra of ^^0(7, p) 
reaction for the rj' meson formation as functions of the 
excited energy which are defined as, 

Eex = m,y - Bri' + [Sp{jp) ~ S'p(p3/2)], (32) 

where B^i is the rj' binding energy and Sp the proton 
separation energy. The rj' production threshold energy 
Ea is indicated in the figure by the vertical dashed line. 
The real part of the optical potential is evaluated for case 
(a) in Eq. (|19|l which has around —150 MeV depth at nu- 
clear center as can be seen in Fig. QJa). The spectra are 
calculated for (upper panel) Wq = — 5 MeV and (lower 
panel) Wo — —20 MeV cases where the imaginary poten- 
tial is written as W[r) = W^pir)/ p^. We also show the 
dominant subcomponents of the spectra in Fig.|Sl As we 
can see from these figures, we can expect to observe the 
peak structures in the spectra due to the formation of 
the ry'-mesic nucleus even in the case with the relatively 
large imaginary potential, and we can expect to deduce 
the magnitude of the rj' mass shift at finite nuclear den- 
sity from the observed spectra. The evaluated imaginary 
part of the ?7'-nucleus potential is small enough and the 
resonance peaks are expected to be clearly separated each 
other. Hereafter, we show the calculated spectra of the 
Ty'-mesic nuclei formation only for Wq = —20 MeV cases. 

In Figs, ini and [7| we show the (7,p) spectra for wider 
energy region including both 77 and 77' production to see 
the anomaly effects in the whole (7,p) spectra simulta- 
neously. To provide realistic spectra, we include the con- 
tributions due to w meson production processes, which 
exist between rj and nj' contributions. As the cj-nucleus 
interaction, we apply the results obtained by two theory 
groups reported in Refs. ^3 and We should men- 
tion here that the w-nucleus interaction at threshold in 
Ref. is attractive, while that in Ref. ^3 is repulsive. 
Hence, we use both interactions to evaluate the contri- 
bution of oj meson. 

In Fig. we show the spectra including the contri- 
bution due to the formation of the uj mesic nuclei with 



the attractive potential reported in Ref. |40j . We can see 
the significant mass reduction of the 77' meson in nuclear 
medium in Fig. El^a), while in Fig. IHIb) we only see the 
simple quasi-free 77' formation spectrum without Ua{^) 
anomaly effect. In the case with the density dependent 
gz) in Fig. EJc), we can see the mass reduction of both 
the 77 and rj' mesons. Especially, even below the lo pro- 
duction threshold, we can see the finite contributions of 
the 77' mesic nuclei formation to the spectrum as the con- 
sequences of the large 77' mass reduction in the case with 
the density dependent go defined as (c) in Eq. ((T^ . For 
the 77 mesic nuclei production, it is hardly to see the dif- 
ferences between the cases with (a) the constant go and 
(c) the density dependence go, due to the large momen- 
tum transfer for the 77 mesic nuclei formation with the 
incident photon energy E^ = 2.7 GeV considered in this 
paper. At this energy, the contributions from higher par- 
tial waves of quasi-free 77 meson production are so large 
that those from the bound 77 states hardly affect the to- 
tal spectrum. Furthermore, we should make a comment 
on the 77 mesic nuclei that, in order to see the anomaly 
effect clearly, we have not considered any contributions 
from resonances like A''* (1535) in this paper in the esti- 
mation of the real part of the 77-nucleus optical poten- 
tial, whose contribution is known to be important on 
the 77-nucleon system. The detailed discussions of the 77- 
mesic nuclei with the A"* (1535) resonance are reported 
in Refs. [s^ [g^. The calculated contribution of the u 
channel in Fig.Ois essentially the same as that reported 
in Ref. |67j|. 



Similarly, in Fig. [7| we show the spectra with the uj 
mesic nuclei with the repulsive potential predicted in 
Ref. li^l. In this case, the contribution of the quasi- 
free LO formation locates the same energy region with the 
bound 77' in the cases (a) and, especially (c) in Fig.[71 In 
Fig. [TJ^b) , we can only see simple three peaks correspond- 
ing the quasi-free 7;, lo and 77' mesons in the spectrum. 



As for the background, which is very important to dis- 
cuss experimental feasibility, is evaluated by using the ex- 
perimental data taken by LEPS collaboration at SPring-8 
recently j?^ . That was a test experiment in the prepara- 
tion stage for the observation of the ui mesic nuclei by the 
(7,p) reaction, which used the same kinematics proposed 
in this paper and observed the background proton emis- 
sion rate from Carbon target including energy region for 
the 77' meson production |73l |. We can roughly estimate 
the order of magnitude of the background proton cross 
section as to be 10 - 100 [nb/sr MeV] in the 77' formation 
region. Thus, we estimate the signal over noise ratio is 
around S/N ^ 1/10. We think the absolute magnitude of 
the calculated formation cross section is reasonably large 
and the spectra are expected to be observed in future 
experiments at SPring-8 j73i |. 
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IV. CONCLUSIONS 

We have investigated the possible observation of the ef- 
fective restoration of the ^/^(l) anomaly at finite density 
by formations of the rj- and ry'-mesic nuclei. Especially, 
since the heavy ry' meson mass is considered to be gen- 
erated from the J7yi(l) anomaly effect, we can expect to 
observe clearly the anomaly effect in medium in the 77' 
channel. 

In order to evaluate more quantitative magnitude of 
the meson mass changes due to the medium effects and 
to obtain deeper insights on J7a(1) anomaly effects than 
Ref. Is^l, we adopted the Nambu-Jona-Lasinio (NJL) 
model with the Kobayashi-Maskawa-'t Hooft (KMT) in- 
teraction, and calculated the meson masses in finite den- 
sity with three different coupling strengths of the KMT 
interaction. We obtained significant 77' mass changes even 
at normal nuclear density due to the effective restora- 
tion of the t^A(l) anomaly. In the finite density calcu- 
lations, we have substituted the SU(2) symmetric con- 
stituent quark matter for the symmetric nuclear matter. 
Although the NJL model doesn't confine quarks, it is 
a good effective model which can describe the low en- 
ergy QCD phenomena, i.e., the properties of light pseu- 
doscalar mesons rather well after fixing the model param- 
eters appropriately. Even the NJL model has some short- 
comings, we think it is meaningful to apply the model 
in this exploratory level to clarify the physics of 7/(958) 
mesic nuclei and Ua{^) anomaly at finite density. 

To investigate the experimental feasibility, we calcu- 
lated formation cross sections of the 77- and yy'-mesic nu- 
clei with (7,p) reaction. We found that the calculated 
cross section has reasonably large magnitude, and the 



(7,p) reaction with GeV photon beam, which can be pro- 
vided in existing facilities like SPring-8, is an appropriate 
reaction for our purpose. We conclude that we can expect 
to observe the 77' mass reduction clearly in this reaction, 
and to obtain new information on the ?7a(1) anomaly at 
finite density. 

The present evaluation is the first theoretical results 
for the formation reaction of the 77- and Ty'-mesic nuclei 
based on the NJL model results to know the behavior 
of C/a(1) anomaly in the medium. We believe that the 
present theoretical results is much important to stimulate 
both theoretical and experimental activities to study the 
C/a(1) anomaly at finite density and to obtain the deeper 
insights of QCD symmetry breaking pattern and the me- 
son mass spectrum. 
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